Wave Functions of Bosonic Symmetry Protected Topological Phases 
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We study the structure of the ground state wave functions of bosonic Symmetry Protected Topo- 
logical (SPT) insulators in 3 space dimensions. We demonstrate that the differences with conven- 
tional insulators are captured simply in a dual vortex description. As an example we show that a 
previously studied bosonic topological insulator with both global U(l) and time-reversal symme- 
try can be described by a rather simple wave function written in terms of dual "vortex ribbons" . 
The wave function is a superposition of all the vortex ribbon configurations of the boson, and a 
factor ( — 1) is associated with each self-linking of the vortex ribbons. This wave function can be 
conveniently derived using an effective field theory of the SPT in the strong coupling limit, and 
it naturally explains all the phenomena of this SPT phase discussed previously. The ground state 
structure for other 3d bosonic SPT phases are also discussed similarly in terms of vortex loop gas 
wave functions. We show that our methods reproduce known results on the ground state structure 
of some 2d SPT phases. 
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I. INTRODUCTION 

The disordered ground states of strongly interacting 
quantum many-body systems can have much richer struc- 
tures compared with classical disordered states. The 
quantum richness of a system is encoded in the entan- 
glement of its ground state wave function, and without 
assuming any symmetry of the Hamiltonian, the ground 
state wave function of a quantum many-body state can 
have long range entanglement, which implies that the 
system has a "topological order" . In the last few years, 
motivated by the discovery of free fcrmion topological in- 
sulators protected by time-reversal symmetry^, it was 
realized that short range entangled state can still be fun- 
damentally distinct from trivial product states, as long as 
the system preserves certain global symmetry G. These 
nontrivial quantum disordered phases with short range 
entanglement are called "symmetry protected topologi- 
cal" (SPT) phases. They are separated from the trivial 
product state through sharp quantum phase transitions 
in the bulk, either continuous or first order. In space di- 
mension d = 1 the Haldane spin chain provides an old and 
nice example of an SPT phased. It has a bulk gap and 
no fractional excitations but nevertheless has dangling 
symmetry protected spin- 1/2 moments at the edge^r— . 
The Haldane chain thus provides an early example of an 
interacting topological insulator. 

In this paper we are mainly concerned with three di- 
mensional symmetry protected topological insulators of 
bosons/spin systems. A formal mathematical classifica- 
tion^ of SPT phases based on group cohomology allows 
a number of such phases to exist (depending on the global 
symmetry) but sheds little light on the physical proper- 
ties. The latter have been discussed recently in Ref. [id]. 
A characteristic feature of all SPT phases is the pres- 
ence of non-trivial surface states at the interface with 
a trivial insulator. Indeed though the bulk is gapped 



and has no fractional excitations or topological order, 
such an interface cannot be in a trivial insulating state. 
Ref. [HI described the effective surface theory for a num- 
ber of three dimensional bosonic topological insulators 
and determined the structure of the allowed non-trivial 
phases. These either spontaneously break the defining 
global symmetry or if gapped have surface topological or- 
der. Exotic symmetry preserving gapless states were also 
shown to be possible. A key feature is that the surface 
effective field theory realizes symmetry in a manner not 
possible in strictly two dimensional systems. Bulk topo- 
logical field theories and effective field theory descriptions 
have also been provided. 

In this work we will study the structure of the ground 
state wave function of various such 3d bosonic SPT in- 
sulating phases with global U(l) and time reversal {Zj) 
symmetries. The differences with conventional Mott in- 
sulators are conveniently captured in a dual description 
in terms of closed vortex loops. In Mott insulating phases 
(conventional or topological) the vortex loops have pro- 
liferated and the ground state wavefunction can be de- 
scribed as a vortex loop gas (see Sec. [TTJ) . We show that 
when compared with the conventional insulator this vor- 
tex loop gas has extra phase factors depending on the 
topology of the vortex loop configuration. We demon- 
strate that these wave functions simply capture all the 
major phenomena of the SPT phases, both in the bulk 
and at the boundary, that were discussed in Ref. [TJ. As 
a key example, in Sec. Mil we discuss a non-trivial SPT 
phase with symmetry either direct (U(l) x Z T ) (as is 
appropriate for a spin model realization of an interact- 
ing boson system) or semidircct (U(l) x Z£) product. 
Here the vortex lines should be viewed as ribbons with 
a non-zero thickness and there is a phase —1 associated 
with each self-linking of a vortex ribbon. We briefly also 
discuss a different SPT phase that occurs for U(l) x Z T 
where each vortex loop can be viewed as a Id Haldane 
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FIG. 1: The wave function of the 3d bosonic SPT discussed in 
this paper is a superposition of all the configurations of vortex 
ribbons with factor (—1) associated with each self- linking. 



spin chain. These results are obtained by analysing both 
the sigma model effective field theory and the topological 
"BF" effective field theories proposed in Ref. [TJ for these 
phases. 

In 2d a result with a similar flavor has been derived 
by Levin and Gv^ for an SPT phase with Ising, i.e Z2, 
symmetry in terms of a domain wall loop gas with phase 
factors. In Sec. IIVI we reproduce this result using our 
methods. We also discuss the ground state wavefunc- 
tion structure of the 2d boson topological insulator with 
U(l) x Z2 symmetry. We use this to obtain a dual vor- 
tex description of this state, and show that the physics 
is correctly captured. 



II. WAVE FUNCTION OF TRIVIAL 3D BOSE 
MOTT INSULATOR 

Let us start with briefly reviewing the trivial Mott in- 
sulating phase of bosons. This is conveniently modeled 
by a quantum disordered phase of interacting U(l) ro- 
tors on a 3d lattice, which is described by the Hamilto- 
nian H = J2<i j> —tcos(9i — 9j) + U{hi) 2 . The boson 
creation operator hi = e lSi and rij is the corresponding 
U(l) charge at site i. 9i and rii are canonically con- 
jugate. The quantum disordered phase of the rotors is 
equivalent to the familiar Mott insulator phase and oc- 
curs when t/U C 1. In the strong coupling limit t — > 0, 
the ground state wave function is a trivial direct product 
state: 

i*>=ni** = °>~n r de ^- w 

Jo 

The wave function of the quantum disordered phase with 
finite but small t/U can be derived through perturbation 
on wave function Eq. [TJ For our purposes it is useful to 
consider a simple approximate form of the wave function 



that captures the physics of the Mott phase : 

|¥>~ r\\d9eMY, -#008(0; -^n^), (2) 
Jo <i,i> 1 

where K ~ t/U <C 1. This wave function is a superposi- 
tion of configurations of 9i with a weight that is the same 
as the Boltzman weight of the 3d classical rotor model. 
The standard duality formalism of the 3d classical ro- 
tor model leads to the dual representation of this wave 
function: 

I*)- J DAj2 e M~ J d 3 x^(V x A) 2 +i2nA- J] 
f 

x \l(x),J(x)}, (3) 

Vector field J takes only integer values on the dual lattice, 
and it represents the vortex loop in the phase 9. In order 
to guarantee the gauge invariance of A, J must have no 
source in the bulk: V • J = 0. The vortex loop J can only 
end at the boundary, which corresponds to a 2d vortex. 
The U(l) gauge field A induces long range interactions 
between vortex loops with coupling strength K . In the 
limit K — > 0, i.e. the strong coupling limit of the original 
rotor, the wave function Eq. [3] for quantum disordered 
lattice bosons becomes a equal weight superposition of 
all vortex loop configurations, with a weak long range 
interaction. 

Quite generally the Mott insulating phase is obtained 
when the vortex loops have proliferated. Consequently 
the ground state wave function can be described as a loop 
gas of oriented interacting vortex loops. The discussion 
above provides a derivation of this loop gas wave function 
starting from a simple but approximate microsopic boson 
wave function. A crucial point about the structure of the 
loop gas wave function for the trivial Mott insulator is 
that it has positive weight for all loop configurations. 



III. WAVE FUNCTION OF 3D BOSONIC SPT 
PHASES 

A 3d SPT phase with U(l) symmetry is also a quantum 
disordered phase of rotor 9i, thus it is expected that its 
wave function is still a superposition of vortex loop con- 
figurations. However, more physics needs to be added 
to the vortex loops in order to capture the novel physics 
of the SPT phase. One of the central results of this pa- 
per is to determine the structure of this vortex loop gas 
wave function for the 3d SPT phases with U(l) and time- 
reversal symmetry discussed in Ref. [Til We first focus 
on one example which occurs for both U(l) x and for 
[7(1) X Z% '• We show that the ground state is described 
by a superposition of vortex loop configurations \C V ), but 
each vortex loop should be viewed as a "ribbon" rather 
than a line, and a self-linking of this ribbon contributes 
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exactly factor (—1) (Fig. [T]): 

m~"£(-l) N ^o{C v }\C v ), (4) 

c„ 

where N t is the number of self-linkings. Here V'o [Cv] is 
the weight of that vortex configuration in a trivial Mott 
insulator. The self-linking of a vortex ribbon is the link- 
ing number between the loops defined by the 2 ends of 
the ribbon. 

This wave function Eq. |4] explains the key phenomena 
of the 3d SPT phase discussed in Ref. [3 In Ref. GJ, 
using a 2+ Id boundary field theory, it was proved that 
the vortex of the U(l) rotor at the boundary of this SPT 
is a fermion. A vortex at the boundary is the end (source) 
of the vortex ribbon in the bulk, and as was discussed 
in Ref. [TH, exchanging the ends of ribbons is equivalent 
to twisting one of the ribbons by 2n, which according 
to Eq. @] should contribute factor (—1). Thus the bulk 
wave function Eq. @] already implies that the vortex at 
the boundary must be a fermion. 

The wave function Eq. 2] can be derived either using 
a bulk non-linear sigma model effective field theory for 
the SPT phase or using a bulk topological "BF" field 
theory. We will present both these derivations below. 
In order to describe the 3d SPT with either U(l) X 
or U(l) x Z 2 symmetry, Ref. [3 proposed the following 
nonlinear Sigma model that involves a five component 
unit vector n — (n , • • • , n 5 ), with a topological O— term 
at 9 = 2tt: 

S = fd 3 xdT -(d»n? + ^-e abcde n a d x n b d v n c d z n d d T n e ,(5) 

J g "4 

where f^4 is the volume of a four dimensional sphere with 
unit radius. Eq. [5] has an enlarged SO(5) symmetry, 
but later we will reduce this symmetry down to physi- 
cal 17(1) x Zj or C/(l) x Zj. ' 

In 3+ld, an order-disorder phase transition occurs 
while tuning g. We will focus on the quantum disordered 
phase with strong coupling g. Since O = 2-7T in Eq. [5l its 
quantum disordered phase has the same bulk spectrum 
as the case with = 0. Thus coupling constant g flows 
to infinity in the quantum disordered phase and this is 
the limit we will focus on in this paper. 

The physical meaning of the 0— term in a NLSM is 
usually interpreted as a factor exp(iO) attached to every 
instanton event in the space-time. Then this interpreta- 
tion would lead to the conclusion that = 27r is equiva- 
lent to = 0. However, this interpretation is very much 
incomplete, because it only tells us that theories with 
= 2n and have the same partition function when 
the system is defined on a compact manifold. These two 
theories actually have very different ground state wave 
functions. In order to expose the wave function, we need 
to keep an open boundary of time. In this case the wave 
function can be derived using the following path integral: 



~ / Dn(x,T)exp(-S)r l ^ =+0C =R',R T= ^ 00 =a- (6) 

The ground state wave function |\E') can then be obtained 
straightforwardly in the strong coupling limit g — > +oo: 

|*) - J Dn{x)W[n]\n(x)) 

Win] — e ixt f d3x So du e abcd< ! n a d3 : n b dyn c dzn d d u n' /^x 

Here ri(x, u) is an extension of the real space config- 
uration n(x) that satisfies H(x, 0) = (0,0,0,0,1), and 
n(x, 1) = n(x). Eq. [7]is a superposition of all the config- 
urations of the 0(5) vector field n(x), with a weight that 
is proportional to the real space Wess-Zumino-Witten 
(WZW) term W[ft] at level- 1. Thus the ground state 
wave function of Eq. [5] with = 2tt is fundamentally 
different from the case with = 0. 

Now let us reduce the artificial SO (5) symmetry of 
Eq.[5]to U(l) and We decompose the five component 
vector n as n — (sm(a)(j>, cos(a)0o), where is a unit 
four component vector, and 0o — ±1 is an Ising order 
parameter. We further define two bosonic rotor operators 
bi ~ + i(f>2, 62 ~ 4>3 + i4>4- Under the U(l) and Zj, 
we take these variables to transform as 

Zl : 61,63 -> 61,62 (C/(l)xZ 2 T ), 

6i,6 2 ->-6l,-&5 (U(l)xZl), 

U(l) : &i->e ie 6i, b 2 ^e l % 2 . (8) 

We assume the system favors 4> over tfio. If the time- 
reversal symmetry is preserved, namely (0o) = 0, the 
WZW term in the wave function Eq.[7]reduces to a theta 
term for the 4-component unit vector ip in 3 + dimen- 
sions. Thus we get the following wave function: 

I*) ~ J D$(x) 

x exp(J ^x^eabcde^rd^d^d^Mx)) 

= Jd$(x)(-1) n °\$(x)). (9) 

This wave function is a superposition of all configurations 
of (j}( x ) in real space, with a 9— term defined in 3d real 
space, at precisely 9 = tt. N s is the Skyrmion number of 
the four component vector cf>, since ^[S* 3 ] = Z. The fact 
9 = it in the wave function is protected by time-reversal 
symmetry Z 2 ■ If this Z 2 symmetry is broken, 9 in this 
wave function will be tuned away from tt. 

For our purposes we need to introduce anisotropies 
that reduce the symmetry from 0(4) to J7(l) x U(l). 
Then the 9 term (at 9 = tt) implies that there is a phase 
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FIG. 2: (a). When the symmetry is (7(1) x (7(1) the bulk wave 
function is a superposition of two flavors of vortex loops with 
factor (—1) attached to each linking, (b—f), braiding between 
two flavors of vortices at the boundary effectively creates one 
extra linking to the bulk vortex loops, which according to 
the bulk wave function would contributes factor (—1). This 
implies that the two flavors of vortices at the boundary have 
mutual semion statistics. 



factor —1 each time the vortex loops of the two boson 
species linkil. Thus this two species boson Mott insu- 
lator has a wave function which is a superposition of all 
vortex loops of the two species with a crucial factor of 
(— 1) L where L is the total number of linked opposite 
species vortex loops. In contrast for the trivial Mott in- 
sulator of the two boson species system, the weight for all 
vortex loop configurations can be taken to be positive. 

It is implicit in the discussion in terms of a four com- 
ponent unit vector <f> that classical configurations of the 
b\ 2 fields are always such that &1.2 cannot simultane- 
ously vanish. As the amplitude of either of these fields 
vanishes in their vortex core this implies that the vortex 
loops of the two species cannot intersect. Thus a config- 
uration with a linking of the two vortex loops cannot be 
deformed to one without a linking. 

This bulk wave function Eq. [9] also implies that at the 
2d boundary, the vortex of b\ and vortex of &2 (sources 
of vortex loops) have a mutual semion statistics, because 
when one flavor of vortex encircles another flavor through 
a full circle, the bulk vortex loops effectively acquire one 
extra linking (Fig. [5]) , which according to the bulk wave 
function would contribute factor (—1). 

Let us now provide an alternate derivation of this result 
using the bulk topological BF theory for the SPT phase 



also proposed in Ref. [l4|. This theory takes the form 

2tt£ 3D = ]T e^Bl v d x ai + £ d,ald x a J a 
1 1, J 

. (10) 

Here B 1 is a rank-2 antisymmetric tensor that is re- 
lated to the current of boson of species 7=1,2 through 
ft = ■h e ^^<yd u B{ lJ . a 1 ^ is a 1-form gauge field which de- 
scribes the vortices of the bosons. Specifically the mag- 
netic field lines of a 1 are identified with the vortex lines 
of the boson of species I. For the SPT state of interest 
the K matrix is simply a x . The parameter = ir (not to 
be confused with the theta parameter in the sigma model 
description) . The crucial difference with the trivial Mott 
insulator is the second 8 term. To get the ground state 
wavefunction we again evaluate the Euclidean path inte- 
gral with open temporal boundary conditions. Using the 
well known fact that the term is the derivative of a 
Chern-Simons term we end up with the following ground 
state wave functional: 

* [alBjg] ~ e <i^/*W*"»J<W^ [a?,fl#] (11) 

Here ipo is the wave functional for the trivial Mott insu- 
lator. The wave functional for the SPT insulator is thus 
modified by a phase factor given by a 3 + dimensional 
Chern-Simons term. As is well known the Chern-Simons 
term is related to a counting of the total linking num- 
ber of the configuration of the magnetic flux lines of the 
gauge fields. Specializing to the case at hand we see that 
in the presence of a 27r flux line of a\, there is a phase 
factor = 7r whenever a 2-7T flux line of ai links with 
it. Thus we reproduce the result that there is a phase of 
7T associated with each linking of opposite species vortex 
lines. 

Finally if the U(l) xU(l) symmetry is broken down to 
diagonal U(l), then the vortex loops of the two species 
will be confined to each other. The resulting common 
vortex loop of the rotor b ~ 61 ~ 62 becomes a ribbon, 
whose two edges are the vortex loops of b\ and 62- Fur- 
ther for simplicity we assume that there is an energetic 
constraint at short distances that prevent two vortex lines 
of the same species from approaching each other. In par- 
ticular we assume that the binding length scale of the 
opposite species vortex loops is smaller than the allowed 
separation between same species vortex loops. Then 
the vortex ribbons cannot intersect each other. Such 
a "hard-core" constraint on the short distance physics 
should not affect the universal long distance behavior of 
the wavefunction^. Note that the binding of the two 
species of vortex loops gives a physical implementation 
of the mathematical concept of 'framing' used to describe 
the topology of knots. The linking between the two fla- 
vors of vortex loops becomes a self-linking of the ribbon. 
Thus wave function Eq.[5]reduces to wave function Eq.2) 
As we mentioned before, this bulk wave function Eq. [4] 
implies that the end point of a vortex ribbon at a 2d 
boundary is a fermion. Similarly bulk external sources 
for vortex ribbons will also be fermions. 
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All these results concur with the boundary theory dis- 
cussed in Ref. [3 There a boundary field theory for the 
SPT is derived, which is a 2+ld NLSM with four compo- 
nent vector <f), and there is a 2+ld space-time O— term at 
precisely O = 7r. This space-time O— term implies that 
the vortex of 62 carries 1/2— charge of 61, and vice versa. 
Thus vortices of b\ and 62 have a mutual semion statis- 
tics. When the symmetry is broken down to one single 
U(l), the 2d vortex at the boundary becomes a bound 
state of the two flavors of vortices: thus eventually this 
bound vortex becomes a fermion. 

Ref. [I!] also described a different interesting SPT phase 
with U(l) x Z2 symmetry. There the surface theory is 
such that the surface vortex carries a Kramers doublet 
in its core. A bulk effective field theory of this phase 
is also obtained^ by starting with the 0(5) non-linear 
sigma model (Eq.[5]) with anisotropies but with a different 
realization of symmetry from the one described above. 
For instance, we can decompose five component vector n 
in a different way: n — (Re[&], Im[6], N x , N y , N z ), where 
6 is a rotor field that transforms under Z^: b — > —b*. 
N is a three component vector that changes sign under 
Z2 but is uncharged under the global U(l). The 3+ Id 
0— term in Eq. [5] implies that the vortex loop of b is in 
a Id Haldane phase of vector N. We may now break the 
symmetry down to just U(l) x Z^. Then a vortex at the 
boundary must carry a Kramer's doublet, because it is 
effectively the edge of the Id Haldane-like phase (with 
Z2 symmetry) along the vortex loop. This is the "Phase 
1 SPT" with symmetry U(l) x Z% discussed in RefA 

Let us now consider the ground state wave function 
for this phase which in terms of the 5-component vector 
is still given by Eqn. [7j Now the interpretation of the 
WZW term is different. As is familiar from discussions of 
deconfined quantum criticality in terms of sigma models 
with WZW termsi£~— in 2 + 1 dimensions, with these 
symmetries an external source for a vortex line carries 
spin-1/2 of the 0(3) rotation that acts on the N vector. 
If the 0(3) symmetry is broken down to Zj, then the 
vortex source still has a Kramers doublet. This of course 
is completely consistent with the picture that each vortex 
may be viewed as a Haldane chain. Thus the ground state 
wave function in this case can be viewed as a vortex loop 
gas of Haldane chains. 



IV. GROUND STATE WAVE FUNCTION OF 2D 
SPT PHASES 

A. 2d SPT phase with Z2 symmetry 

Let us now switch gears to SPT phases in 2d. We be- 
gin by making contact the work of Levin and Gn-^- on 
the Ising SPT phase. The simplest SPT phase in 2d has 
a Z2 global symmetry. In Ref. |l5|, a lattice model for 
this phase has been discussed. The ground state wave 
function was argued to be a sum over all domain wall 




FIG. 3: (a) Skyrmion of 0(3) vector ft in 2d space. (6) If the 
SO(3) symmetry is broken down to Z2, the Skyrmion becomes 
a domain wall of Z2 order parameter n z . 

configurations with a factor (— 1) d where Nd is the to- 
tal number of domain wall loops. We now show how to 
reproduce this results within the methods of this paper. 

In 2 + 1 space-time dimensions many SPT phases are 
conveniently described by starting with an effective non- 
linear sigma model field theory in terms of a four com- 
ponent unit vector <p. The field theory action reads 

S= [ d 2 xdr -(dj) 2 + ^eabed^p^d^dv^dp^) 
J g 127H 

The crucial ingredient is the term for the 4-component 
unit vector at = 27r. This action Eq. [T^Jhas an SO(4) 
symmetry, and this SO (4) symmetry contains a subgroup 
Zi symmetry <j) — > — cf>. Eventually we will break the ar- 
tificial SO(4) symmetry of Eq. [T2l down to this Z2 sub- 
group. 

In the paramagnetic phase, i.e. in the limit g — > +00, 
the bulk ground state wave function is 

f - i2n 
~ y D<j>(x) exp{^ x 

Jd 2 xj due abcd e^ p rdp4> b dv<p c d p ^ d } |0(x))(13) 

which is a superposition of all configurations of 4>{x) with 
a real space WZW weight. 

Now let us decompose the four component vector into 
4> = (cos(a)0o, sin(a)n), where (f>Q = ±1 is an Ising 
order parameter, and n is a unit three component vector: 
(ft) 2 = 1. We also break the SO (4) symmetry down to 
Z 2 x SO (3) symmetry: 

Z2 ■ 4>o -> ~4>o, n -> -n, 
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S0(3) : Rotation of ft. (14) 

Under this symmetry reduction, if the system energeti- 
cally favors vector ft over <j) (favors a — 7r/2), the wave 
function Eq. [T3] reduces to 

*> ~ J Dn{x) ex P{^ J d 2x £abc£ t ivn a d t _ l n b d„n c } \n(x)} 
= ^(-l)^|n(a:)), (15) 

where N s is the number of Skyrmions of 0(3) vector ft 
in the 2d space. As long as we keep the Z 2 symmetry 
<t> — > —<f>, the expectation value of 4>q is zero, then each 
Skyrmion will contribute a phase factor exactly (—1). 

Now let us further break Z 2 X SO (3) symmetry down 
to Z 2 x S0(2): 

Z 2 ■ 4>o -> -0o, n -> ~n, 
SO(2) : Rotation of n x , n v . (16) 

Also we assume that the system favors n z over n x and n v , 
then each skyrmion becomes a domain wall of Z 2 order 
parameter n z (x) (Fig. [3]) . Now the wave function Eq. [15] 
reduces to a superposition of configurations of Z 2 order 
parameter n z (x): 

(17) 

where is the number of closed domain wall loops of 
n z {x). Eventually we can also break the residual SO(2) 
symmetry, and the wave function Eq. [T7]is unchanged. 

The wave function Eq. [17] is exactly the one derived 
from the lattice model of 2d SPT phase with Z 2 symme- 
try^. In the appendix we will also demonstrate that the 
effective field theory Eq . [T2l implies that, after coupling 
this SPT phase to a dynamical Z 2 gauge field, the 7r— flux 
of this Z 2 gauge field has a semion statistics, which is 
consistent with the result in Ref. [l5|. 

With a full SO(4) symmetry, the edge states of Eq.[T2l 
with precisely O = 2ir is the nonchiral SU(2)! conformal 
field theory (or equivalently as an SO (4) non- linear sigma 
model with a level- 1 WZW term). Since the original 
SO (4) symmetry is reduced to its Z 2 subgroup <j> — > — d>, 
we have to argue that the edge state of Eq. [T2] survives 
under this symmetry reduction. Because the Z 2 sym- 
metry acts on all four components of <f>, in the bound- 
ary WZW model, terms allowed by the Z 2 symmetry are 
J2i j gij4>i4>j (i, J — 0, 1, 2, 3). If these terms are relevant, 
it leads to spontaneous Z 2 symmetry breaking and two 
fold degeneracy at the boundary. Thus the edge state 
cannot be completely trivial (gapped and nondegenerate) 
as long as the Z 2 symmetry is preserved. 



B. 2d SPT phase with (7(1) x Z% symmetry 

Finally we consider the 2d bosonic topological insula- 
tor which occurs when the global symmetry is U(l) x Z 2 . 
This may be described by starting again with the same 
4-component non-linear sigma model but with the follow- 
ing implementation of the physical symmetry. We write 
<f>2— i4>3 = b and let b have charge- 1 under the global U(l). 
Under Z 2 , we demand b — > b, (j>o — > —<fio, <fii — > —<t>i- As 
before we again assume first an anisotropy that prefers ft 
over 0o so that the ground state wave function is given 
by Eqn. [151 Now we introduce further anisotropy to 
reduce to the desired U(l) x Z 2 . The defects of the 
charged field b are of course point vortices. In the core of 
these vortices the amplitude of b is suppressed and the cf) 
points entirely in the <p\ direction. There are two differ- 
ent vortices - known as merons - depending on whether 
in the core <pi = ±1. Each meron may be viewed as 
half a skyrmion and has N s = sgn(0i)l/2. Thus the 
ground state wave function is then a sum over all pos- 
sible configurations of the two kinds of meron vortices 
with phase factors e ±z ^ for the two kinds of vortex. Let 
n v ± be the vortex number of either species at site i in a 
lattice description. Then we require that the total vor- 
ticity N v = n v+ + n„_ = 0. Then the phase factor in 
the wave function = Thus 

there is a relative phase of — 1 associated with — vortices 
compared to + vortices. 

We now argue that this structure of the wave function 
matches what is known about the 2d bosonic topological 
insulator. Consider a dual description of such an insu- 
lator. From our arguments above there are two kinds of 
vortex fields <& v ± corresponding to the two meron vor- 
tices. The dual vortex theory will have a Lagrangian 

C d =^2\{d„- ia^ vs \ 2 + ... + | (e M „AcU A ) 2 (18) 

Here a M is the usual fluctuating non-compact f (1) gauge 
field whose flux density is the original boson number. 
There must in addition be terms where the meron cores 
tunnel into each other: 

A + fc.c) (19) 

We begin by ignoring these and we will reinstate them 
later. Under time reversal a vortex must go to an anti 
vortex (as the boson phase is odd) and the meron cores 
flip into each other. Thus under Z 2 

$„+ -> ±$1_ (20) 

In the trivial insulator all vortex configurations con- 
tribute with the same sign and we must choose $.„ + — >• 
+$^_ . But for the topological insulator there is a relative 
— sign between the two vortex species. Thus we must 
choose $„ + — > — 3>J_. Condensing vortices that trans- 
form in this manner will give us the boson topological 
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insulator. Now let us include the meron core tunneling 
term. Then the two vortex species mix with each other 
so that we identify a single vortex $ w = ~ Its 
transformation under time reversal is & v — > ±$J where 
the + sign describes the trivial insulator and the — the 
topological insulator. This is exactly the same transfor- 
mation law for the vortices that is dictated by the edge 
theory analysis of the 2d boson topological insulato r 20 ' 21 . 
Thus the wave function description we developed cap- 
tures the physics of this state, and further gives a bulk 
dual vortex description. 

V. SUMMARY 

In summary, we have demonstrated in this work that, 
although most of the novel phenomena of a SPT phase 
occur at its boundary, its bulk ground state wave func- 
tion is indeed drastically different from a trivial direct 



product disordered phase. This bulk wave function can 
be conveniently derived from the effective field theory of 
the SPT phase. The structure of the ground state wave 
functions in terms of dual vortex configurations derived 
in this work provide a simple physical picture of the phe- 
nomena associated with these SPT phases. 
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Appendix A: Appendix: Dynamical Z-z gauge fields 
in the 2d Ising SPT 

In this appendix we demonstrate that the effective field 
theory Eq. [T^] not only gives us the correct ground state 
wave function (Eq. [T7|) of the 2d Ising SPT phase, after 
coupling the SPT phase to a Z 2 gauge field, the topolog- 
ical 0— term of Eq. [T2]also leads to nontrivial statistics 
of the dynamical 7r— flux (vison) of the Zi gauge field. 

First of all, Eq. [T^Jcan be rewritten as a SU(2) princi- 
ple chiral model, by introducing SU(2) matrix field G — 
4>°a° +i$-a. G has SU(2)-left and SU(2) -right trans- 
formations: G -> vlGV R . Let us "gauge" SU(2)-left 
and SU(2)— right transformations with dynamical U(l) 
gauge fields a^a z and & M er 2 , i.e. replace d^G with dpG + 
ia^a z G + ib^Ga*. According to Ref. I22II23I . after inte- 
grating out matrix field G, gauge fields and b^ both 
acquire a Chern-Simons term: 

d 2 xdr —e^u P a^d v a p - —e^pb^bp. (Al) 

This is because Eq . [T2l also describes a U(l) bosonic SPT 
with Hall conductivity 2. 

A dynamical U(l) gauge field with level— A: has the fol- 
lowing properties: its charged quasiparticle carries gauge 
flux 2n/k, and this quasiparticle has a statistics angle 
ir/k. Thus the Chern-Simons action Eq. I All gives the 
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7T— flux of U(l) gauge field a M and a semion statis- 
tics, with statistics angle +ir/2 and — ir/2 respectively. 
Notice that the two U(l) gauge groups share the same 
Z2 transformation G — > — G, thus we can break the two 
U(l) gauge fields down to one Z2 gauge field, then the 
dynamical 7r— flux of this Z2 gauge field has two differ- 



ent flavors with semionic statistics angle +ir/2 and — it/2 
respectively. 

In Ref. 15. using their lattice model, the authors con- 
cluded that the dynamical 7r— flux of this Z2 gauge field 
has a semion statistics. Here we have derived the same 
result using our field theory Eq. [T2l 



